New indicators for monitoring genetic diversity applied to alpine brown trout populations using whole genome sequence data.

International policy recently adopted commitments to maintain genetic diversity in wild populations to secure their adaptive potential, including metrics to monitor temporal trends in genetic diversity - so-called indicators. A national programme for assessing trends in genetic diversity was recently initiated in Sweden. Relating to this effort, we systematically assess contemporary genome-wide temporal trends (40 years) in wild populations using the newly adopted indicators and whole genome sequencing (WGS). We use pooled and individual WGS data from brown trout (Salmo trutta) in eight alpine lakes in protected areas. Observed temporal trends in diversity metrics (nucleotide diversity, Watterson's Ï´ and heterozygosity) lie within proposed acceptable threshold values for six of the lakes, but with consistently low values in lakes above the tree line and declines observed in these northern-most lakes. Local effective population size is low in all lakes, highlighting the importance of continued protection of interconnected systems to allow genetic connectivity for long-term viability of these populations. Inbreeding (FROH ) spans 10%-30% and is mostly represented by ancient (<1 Mb) runs of homozygosity, with observations of little change in mutational load. We also investigate adaptive dynamics over evolutionarily short time frames (a few generations); identifying putative parallel selection across all lakes within a gene pertaining to skin pigmentation as well as candidates of selection unique to specific lakes and lake systems involved in reproduction and immunity. We demonstrate the utility of WGS for systematic monitoring of natural populations, a priority concern if genetic diversity is to be protected.

Larval transcriptomes reflect the evolutionary history of plant-insect associations.

In this study, we investigated whether patterns of gene expression in larvae feeding on different plants can explain important aspects of the evolution of insect-plant associations, such as phylogenetic conservatism of host use and re-colonization of ancestral hosts that have been lost from the host repertoire. To this end, we performed a phylogenetically informed study comparing the transcriptomes of 4 nymphalid butterfly species in Polygonia and the closely related genus Nymphalis. Larvae were reared on Urtica dioica, Salix spp., and Ribes spp. Plant-specific gene expression was found to be similar across butterfly species, even in the case of host plants that are no longer used by two of the butterfly species. These results suggest that plant-specific transcriptomes can be robust over evolutionary time. We propose that adaptations to particular larval food plants can profitably be understood as an evolved set of modules of co-expressed genes, promoting conservatism in host use and facilitating re-colonization. Moreover, we speculate that the degree of overlap between plant-specific transcriptomes may correlate with the strength of trade-offs between plants as resources and hence to the probability of colonizing hosts and complete host shifts.

Improved chromosome-level genome assembly of the Glanville fritillary butterfly (Melitaea cinxia) integrating Pacific Biosciences long reads and a high-density linkage map.

BACKGROUND: The Glanville fritillary (Melitaea cinxia) butterfly is a model system for metapopulation dynamics research in fragmented landscapes. Here, we provide a chromosome-level assembly of the butterfly's genome produced from Pacific Biosciences sequencing of a pool of males, combined with a linkage map from population crosses. RESULTS: The final assembly size of 484 Mb is an increase of 94 Mb on the previously published genome. Estimation of the completeness of the genome with BUSCO indicates that the genome contains 92-94% of the BUSCO genes in complete and single copies. We predicted 14,810 genes using the MAKER pipeline and manually curated 1,232 of these gene models. CONCLUSIONS: The genome and its annotated gene models are a valuable resource for future comparative genomics, molecular biology, transcriptome, and genetics studies on this species.

Exploring a Pool-seq-only approach for gaining population genomic insights in nonmodel species.

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool-seq data to generate a de novo genome assembly for mining exons, upon which Pool-seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual-based single nucleotide polymorphisms [SNPs]) and from mapping the Pool-seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F ST) between the two introduced populations exceeds that of the naturally sympatric populations (F ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ ( π ¯  ≈ 0.002 and π ¯  ≈ 0.001 when mapping to S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high-quality reference assembly from a divergent species. We conclude that the Pool-seq-only approach can be suitable for detecting and quantifying genome-wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.

Unprecedented reorganization of holocentric chromosomes provides insights into the enigma of lepidopteran chromosome evolution.

Chromosome evolution presents an enigma in the mega-diverse Lepidoptera. Most species exhibit constrained chromosome evolution with nearly identical haploid chromosome counts and chromosome-level gene collinearity among species more than 140 million years divergent. However, a few species possess radically inflated chromosomal counts due to extensive fission and fusion events. To address this enigma of constraint in the face of an exceptional ability to change, we investigated an unprecedented reorganization of the standard lepidopteran chromosome structure in the green-veined white butterfly (Pieris napi). We find that gene content in P. napi has been extensively rearranged in large collinear blocks, which until now have been masked by a haploid chromosome number close to the lepidopteran average. We observe that ancient chromosome ends have been maintained and collinear blocks are enriched for functionally related genes suggesting both a mechanism and a possible role for selection in determining the boundaries of these genome-wide rearrangements.

Why stay in a bad relationship? The effect of local host phenology on a generalist butterfly feeding on a low-ranked host.

BACKGROUND: In plant-feeding insects, the evolutionary retention of polyphagy remains puzzling. A better understanding of the relationship between these organisms and changes in the metabolome of their host plants is likely to suggest functional links between them, and may provide insights into how polyphagy is maintained. RESULTS: We investigated the phenological change of Cynoglossum officinale, and how a generalist butterfly species, Vanessa cardui, responded to this change. We used untargeted metabolite profiling to map plant seasonal changes in both primary and secondary metabolites. We compared these data to differences in larval performance on vegetative plants early and late in the season. We also performed two oviposition preference experiments to test females' ability to choose between plant developmental stages (vegetative and reproductive) early and late in the season. We found clear seasonal changes in plant primary and secondary metabolites that correlated with larval performance. The seasonal change in plant metabolome reflected changes in both nutrition and toxicity and resulted in zero survival in the late period. However, large differences among families in larval ability to feed on C. officinale suggest that there is genetic variation for performance on this host. Moreover, females accepted all plants for oviposition, and were not able to discriminate between plant developmental stages, in spite of the observed overall differences in metabolite profile potentially associated with differences in suitability as larval food. CONCLUSIONS: In V. cardui, migratory behavior, and thus larval feeding times, are not synchronized with plant phenology at the reproductive site. This lack of synchronization, coupled with the observed lack of discriminatory oviposition, obviously has potential fitness costs. However, this "opportunistic" behavior may as well function as a source of potential host plant evolution, promoting for example the acceptance of new plants.

Mechanisms of macroevolution: polyphagous plasticity in butterfly larvae revealed by RNA-Seq.

Transcriptome studies of insect herbivory are still rare, yet studies in model systems have uncovered patterns of transcript regulation that appear to provide insights into how insect herbivores attain polyphagy, such as a general increase in expression breadth and regulation of ribosomal, digestion- and detoxification-related genes. We investigated the potential generality of these emerging patterns, in the Swedish comma, Polygonia c-album, which is a polyphagous, widely-distributed butterfly. Urtica dioica and Ribes uva-crispa are hosts of P. c-album, but Ribes represents a recent evolutionary shift onto a very divergent host. Utilizing the assembled transcriptome for read mapping, we assessed gene expression finding that caterpillar life-history (i.e. 2nd vs. 4th-instar regulation) had a limited influence on gene expression plasticity. In contrast, differential expression in response to host-plant identified genes encoding serine-type endopeptidases, membrane-associated proteins and transporters. Differential regulation of genes involved in nucleic acid binding was also observed suggesting that polyphagy involves large scale transcriptional changes. Additionally, transcripts coding for structural constituents of the cuticle were differentially expressed in caterpillars in response to their diet indicating that the insect cuticle may be a target for plant defence. Our results state that emerging patterns of transcript regulation from model species appear relevant in species when placed in an evolutionary context.

Micro-injection of Lygus salivary gland proteins to simulate feeding damage in alfalfa and cotton flowers.

Alfalfa and cotton flowers were pierced with small glass capillaries of an overall size and shape similar to that of Lygus stylets, and injected with small quantities (6 to 100 nL) of solutions that contained Lygus salivary enzymes. Crude and partially purified protein solutions from Lygus heads and isolated salivary glands showed substantial polygalacturonase (PG) activity, as has been previously reported. Following injection with both crude and partially purified protein solutions, as well as with pure fungal and bacterial PGs, flowers of both alfalfa and cotton exhibited damage similar to that caused by Lygus feeding. Injection with the same volume of a buffer control as well as a buffer control containing BSA at a comparable protein concentration (approximately 6 microg/mL) showed no symptoms. These results are consistent with a previously suggested hypothesis that the extensive tissue damage caused by Lygus feeding is primarily due to the action of the PG enzyme on the host tissue, rather than to mechanical damage caused by the insect stylet. Substantial genotypic variation for a PG inhibiting protein (PGIP) exists in alfalfa and cotton. We, therefore, suggest that breeding and selection for increased native PGIP levels, or transformation with genes encoding PGIP from other plant species, may be of value in obtaining alfalfa and cotton varieties that are more resistant to Lygus feeding damage.